A limit on the detectability of the energy scale of inflation 
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We show that the polarization of the cosmic microwave background can be used to detect gravity 
waves from inflation if the energy scale of inflation is above 3.2 x 10 15 GeV. These gravity waves 
generate polarization patterns with a curl, whereas (to first order in perturbation theory) density 
perturbations do not. The limiting "noise" arises from the second-order generation of curl from 
density perturbations, or rather residuals from its subtraction. We calculate optimal sky coverage 
and detectability limits as a function of detector sensitivity and observing time. 
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Few ideas have had greater impact in cosmology than 
that of inflation |[], |[ Inflation makes four predic- 
tions, three of which provide very good descriptions of 
data: the mean curvature of space is vanishingly close to 
zero, the power spectrum of initial density perturbations 
is nearly scale-invariant, and the perturbations follow a 
Gaussian distribution. As the data have improved sub- 
stantially (e.g., 0, H, El) they have agreed well with in- 
flation, whereas all competing models for explaining the 
large-scale structure in the Universe have been ruled out 
(e.g.,011). 

We must note though that these three predictions are 
all fairly generic pgfl. Further, although existing models 
for the formation of structure have been ruled out, there 
is no proof of inflation's unique ability to lead to our 
Universe. Indeed, alternatives are being invented [ |To| |. 

The fourth (and yet untested) prediction may there- 
fore play a crucial role in distinguishing inflation from 
other possible early Universe scenarios. Inflation in- 
evitably leads to a nearly scale-invariant spectrum of 
gravitational waves, which are tensor perturbations to 
the spatial metric. Detection of this gravitational-wave 
background might allow discrimination between compet- 
ing scenarios (e.g. and different inflationary models 
(e.g. 0). 

The amplitude of the power spectrum of tensor per- 
turbations to the metric is directly proportional to the 
energy scale of inflation. One can use a determination of 
the tensor contribution to CMB temperature anisotropy, 
here parameterized by the quadruple variance, to deter- 
mine this energy scale 

u» 1/4 /™p* = 1-2(Qt) 1/4 = 3.0 x ictV/ 4 (l) 

where r = (Qt)/(Q%): & stands for scalar (density) per- 
turbation and (Q%) - 4 x 10~ n from observations pfj[ . 
Without detection of gravitational waves, the energy 
scale of inflation remains uncertain by at least 12 orders 
of magnitude. Pinning down this energy scale could be 
crucial to understanding how inflation arises in a funda- 
mental theory of physics. 

In Fig. 1 we show the angular power spectrum of CMB 
temperature perturbations contributed by scalar pertur- 
bations and by tensor perturbations with r = 1CP 3 . By 
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FIG. 1: Angular power spectra. Solid lines are for tempera- 
ture anisotropies due to scalar perturbations, C^i and tensor 
perturbations, C^; with r = 10 -3 . Dashed lines are for the 
E modes from scalar perturbations, C§ ; and the B modes 
from tensor perturbations, Cjjj. The dotted lines are for the 
lensing-induced scalar B-modes, C% 1 before (above) and after 
(below) the cleaning that can be done by a perfect experiment. 



determining the total CMB temperature power spectrum 
we can determine or limit the energy scale of inflation, 
based on the presence or absence of extra power at low 
I. The scalar temperature perturbations inevitably limit 
our ability to detect the tensor temperature perturba- 
tions to those cases with r > n im = 0.13. |fl3| [|30|| 

In [|l4| [l5| it was pointed out that tensor perturbations 
result in CMB polarization patterns with a curl, whereas 
scalar perturbations do not. By analogy with clcctromag- 
netism, these modes are called "B-modes" , and the curl- 
free modes are called "E-modes". This was an exciting 
development, because the new signature of tensor pertur- 
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bations could not be confused with scalar perturbations. 
Figure 1 also shows the power spectrum of this B-mode 
with the amplitude it would have if V^ 4 = 6.4 x 10 15 
GeV, corresponding to r = 1CP 3 . It is a very weak 
signal, even at its peak 10 7 times less than the power 
spectrum of the temperature anisotropy. Ignoring any 
contaminating source of curl-mode polarization the de- 
tectability limit is solely a matter of sufficient sensitivity 
to the CMB polarization. For a detector of sensitivity s 
uniformly observing the entire sky for time t: 



Hi: 



10" 



/iKy/sec 



t 



1 year 



(2) 



Equation g [[l6| does not take into account a contam- 
inating source of B mode that arises from the lensing of 
the E mode by density perturbations along lines-of-sight 
between the observer and the last-scattering surface |l7J . 
This scalar contribution to the B-mode power spectrum 
is shown in Fig. 1. As we will see below this contam- 
ination sets the detectability limit at ru m — 6 x 10~ 4 , 
similar to what was found in |Ts|] . 

This lensing contaminant can be cleaned from the 
maps as discussed in |l?J, ^FJ. Here we show that there 
are limits to how well these maps can be cleaned. Our 
main result is that in the no-noise limit the tensor per- 
turbations can be detected only if V» 1/4 > 3.2 x 10 15 GeV. 
Interestingly, this is below the mass scale for the gauge 
coupling-constant convergence at 10 16 GeV that occurs 
in minimal extensions to the standard model of particle 



physics (e.g. (21)). 

The Stokes parameters, /, Q and U are related to the 
unlensed Stokes parameters (denoted with a tilde) by 



1(0) = 1(0 - 



Q(0) = Q(0 + 50) U(0) = U(0- 



The deflection angle, 60, is the tangential gradient of the 
projected gravitational potential, 



( v — v ) 
dr- — ^>(rh, r) 



(4) 



where r is the coordinate distance along our past light 
cone, s denotes the CMB last-scattering surface and n is 
the unit vector in the direction. 

The effect on the B-mode power spectrum is piT 



Cbi = c §1 + J2w!'c §1 , 



(5) 



where 
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oj}(0) lMmMl'6) + J o (l'0)J o (W)}}. (6) 



The functions <J 2 (0) and <Jq(0) depend on the statistical 
properties of the displacement potential <f> and are given 
by 
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and 4>i m is the spherical-harmonic transform of i 

If we have a means of determining </>, and therefore 60, 
we can reconstruct the unlensed maps from the lensed 
maps by use of Eq. ||. This procedure cleans out the 
lensing-induced B-mode. The Cf above can either be 
interpreted as the power spectrum of the lensing potential 
(in the case of uncleaned maps in which case we will call 
it Cf (S) ) or the power spectrum of the lensing potential 
residuals (in the case of lensing-cleaned maps). In the 
uncleaned case [ p2) (in the Limber approximation valid 
at small scales): 
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where r s is the comoving distance to the last-scattering 
surface and A|(fc,r) = fc 3 /(27r) 2 P$(/c, r) where P$(/c,r) 
is the power spectrum of the gravitational potential at 
the time corresponding to coordinate distance r on our 
past light cone. We plot in Fig. 2 

We consider the (j) reconstruction procedure given in 
p0| which exploits the fact that lensing leads to a mode- 
mode coupling with expectation value proportional to <f>. 
Their estimator is a minimum-variance (MV) average 
over pairs of map modes with / ^ V . Since we know the 
statistics of the signal, Cf , we can Wiener filter (WF) 
the MV estimate and further reduce the errors in the 
reconstruction. The error in the WF estimate of 4>i m has 
variance: 
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(10) 



The Wiener-filtering is important in the low signal-to- 
noise regime (Cf /Cf <C 1) where the variance of 
the reconstructed 4> is instead of cf^ MV \ 

Cleaning can greatly reduce the amplitude of Cf as is 
shown in Fig. 2 where the angular power spectrum of 
the residual <p is shown for two different experiments. 
The first is the 'reference' experiment of |2(J which 
makes temperature maps with weight-per-solid angle of 
w = (1/iK — arcmin)~ 2 and Q and U maps each with 
half this weight, all with 7' (full- width-half-maximum) 
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the maps are noisy and the lensing contaminant can not 
be removed accurately enough, too little sky hampers 
the statistical subtraction of the residual scalar B-mode 
power spectrum. The optimal sky coverage is given by 
f°& = l(r 8 (>K) 2 i/s 2 - In Fig. 3 we plot r Um as a func- 
tion of s 2 ft assuming optimal sky coverage (constrained 
to / s ky < 1)- The detectability limit is not strongly 
sensitive to deviations from optimal sky coverage, which 
is fortunate since other considerations can influence sky 
coverage choice as well. Note that s 2 /t = 10 -8 (/iK) 2 
could be achieved by observing for a year with an array 
of 30,000 detectors each with sensitivity 100/iK^/sec. Ten 
times more detectors (or integration time) are required 
to reach s 2 /t — 10~ 9 (/iK) 2 and approach the smallest 
possible Him- We have ignored aliasing of the scalar E- 
mode wich can be important for / s k y ^ 1, as studied in 
Q. Aliasing considerations increase f°£* and rn m . For 
s 2 /t = 10- 7 (^K) 2 and / sky = 1, r lim = 7 x 10~ 4 . There- 
fore, including effects of aliasing would not increase rn m 
by more than 40% over our results as shown in Fig. 3. 



FIG. 2: Angular power spectrum of projected gravitational 
potential, (j> (solid curve) and the power spectrum of the resid- 
uals of different reconstruction procedures: minimum variance 
for "reference" experiment (long dashes) , Wiener filter for the 
reference experiment (short dashes) and Wiener filter for a 
noiseless experiment (dotted line). 



resolution. The second is the no-noise, perfect angular 
resolution case. One can see that high signal-to-noise 
reconstruction of (j> is possible. 

Because the cleaning is not perfect, there is still 
lensing-induced scalar B-mode in the cleaned maps. In 
the no-noise limit the residual scalar B-mode power spec- 
trum, Cg[ F , has been reduced by a factor of 10 from the 
uncleaned amplitude, as shown in Fig. 1. 

Writing the tensor-induced B-mode power spectrum 
as = rC^f the resulting error in r, oy, is given by 

(ii) 

where the Q and U maps cover f s ^ y of the sky with uni- 
form weight per solid angle w and have been smoothed 
by a circular Gaussian beam profile with full-width- 
half- maixmum = \/81n2<7(,. For riim = 3.3oy we find 
for noise-free maps, but with no cleaning done, rn m — 
6 x 10~ 4 . For the 'reference' experiment of we find 
Him = 4 x 10 -4 after cleaning with the MV estimator and 
1.8 x 10~ 4 with the WF estimator. Finally, for noise-free 
maps with WF cleaning r\i m = 6 x 10 -5 . 

All of the above results are for / s k y = 1. Given a 
detector (or array of detectors) with total sensitivity s 
and an amount of observing time t there is an optimal 
amount of sky to cover. Too much sky coverage means 
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FIG. 3: Achievable detectability limit as a function of total 
detector array sensitivity, s and observing time t, assuming 
optimal sky coverage (lower curve) and 1/3 of optimal sky 
coverage (upper curve). Detectability limit assuming three 
times optimal sky coverage overlaps the upper curve when 
optimal sky coverage is less than 1/3. 



Since we are interested in determining a fundamen- 
tal limit to the detectability of tensor perturbations, un- 
deterred by the virtual cost of imaginary experiments, 
we explore one further method of cleaning out the lens- 
ing contaminant. Imagine we reconstruct the projected 
gravitational potential or the density field on our past 
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light cone out to some limiting redshift, zii m jn]- Then 
we could directly calculate the lensing generated by the 
matter at z < zy lm . Subtracting this lensing contami- 
nant from CMB maps leaves only the contribution from 
z > z\im- The effect of this remaining lensing can be 
calculated by altering the lower limit of integration in 
Eq. H from to r(z\ im ). To improve upon the noise-free 
self-cleaned maps result of rn m = 6 x 10~ 5 one needs 
to take Z\i m > 7. We know of no way of obtaining this 
information at such high z. 

Our analysis ignores polarized emission from galactic 
and extragalactic sources. Multi-frequency observations 
can be used to clean out these signals based on their 
distinct spectral shapes. However, even in the no-noise 
limit these cannot be cleaned perfectly since their spec- 
tral shapes are not perfectly well known and vary spa- 
tially (e.g. [^3)). Our estimate of the detectability limit 
should be viewed as a lower limit. 

How likely is it that the energy scale of inflation is high 

1 /4 

enough to be determined? It is conceivable that V* is 
as low as 10 3 GeVQ. Given even odds between 10 3 
and 10 19 GeV the chances are small. On the other hand 
gauge coupling unification gives us a hint that something 
interesting may be occurring at 10 16 GeV. If inflation 
has anything to do with grand unification, or physics at 
higher energy scales, then the chances are good. 

The likelihood of detection may also be addressed by 
better determination n f the scalar spectrum Tf wp a.s- 



Hiimln particular functional forms fnr V(if>) and » scalar 



perturbation spectrum with power-spectral index, n, 
near the scale-invariant value of unity, then we can give 
constraints on the range of possible values of r. Here we 
follow the nomenclature and calculations of p| . Expo- 
nential potentials have r = 5(1— n) and "small-field poly- 
nomial" potentials (V((p) = A 4 [1 - (0/m) P ]) with p = 2 
have r — 10 exp (—50 (1 — n)). Still other classes of mod- 
els (hybrid inflation and the small-field polynomial po- 
tentials withp > 2) leave little or no relationship between 
r and n; r can take on detectable values or vanishingly 
small ones. Perhaps the happiest case is the simplest 
one: polynomial potentials with V((f>) = A 4 ((j)/ /j,) p . With 
p > 2 all have r > 0.1, well above our detectability limit. 

In sum, we have calculated how large the amplitude 
of tensor perturbations must be in order for their effect 
on the CMB to be distinguishable from those of lensing. 
By reconstructing the lensing potential, using the mode- 
mode coupling induced by lensing, one can reduce the 
lensing contamination but not eliminate it. The residual 
lensing signal prevents detection of tensor perturbations 
if V* < (3.2 x 10 15 GeV) 4 . If K is slightly larger then 
an ambitious observational program may succeed in de- 
termining the energy scale of inflation — possibly a key 
step towards finding inflation a comfortable home in a 
fundamental theory of physics. 

We thank A. Albrecht, D. Chung, N. Dalai and J. 
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For example, a much more impressive success would be 
prediction of features in the power spectrum followed by 
their detection. 

The * subscript on V means that it is evaluated when the 
scale factor is about e 50 times smaller than it is at the end 
of inflation, when the relevant perturbations for large- 
scale CMB fluctuations are exiting the Horizon. There is 
a slight dependence on J7a j2(j], here taken to be SI a = 
0.65. Other cosmological parameters assumed throughout 
are h = 0.65, Q b = 0.05, Q m = 0.35 and 8 H = 4.2 x 10~ 5 . 
Here and throughout, the detection limit is set at 3.3 
times the standard deviation. For a normal distribution 
this means that one can make a detection at 95% confi- 
dence for 95% of possible realizations. 
Perhaps this is done with a POTENT-likeQ use of pe- 
culiar velocity measurements, or analysis of the statistical 
distribution of apparent galaxy shapes. 



